Two dimensional XXZ-Ising model on square-hexagon lattice 
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We study a two dimensional XXZ-Ising on square-hexagon (4-6) lattice with spin-1/2. The phase 
diagram of the ground state energy is discussed, shown two different ferrimagnetic states and two 
type of antiferromagnetic states, beside of a ferromagnetic state. To solve this model, it could 
be mapped into the eight-vertex model with union jack interaction term. Imposing exact solution 
condition we find the region where the XXZ-Ising model on 4-6 lattice have exact solutions with one 
free parameter, for symmetric eight- vertex model condition. In this sense we explore the properties 
of the system and analyze the competition of the interaction parameters providing the region where 
it has an exact solution. However the present model does not satisfy the free fermion condition, 
unless for a trivial situation. Even so we are able to discuss their critical points region, when the 
exactly solvable condition is ignored. 



Recently frustrated magnetic systems have been at- 
tracting a lot of attention due to their rich properties. 
Such systems have several phase diagrams displaying a 
number of unusual quantum phases 0, 0|. Frustration 
interaction is exhibit experimentally in inelastic neu- 
tron scattering. Then two dimensional magnetic lattice 
is a challenge for the theoretical investigation. After 
Onsager'sQ solution of the square two dimensional Ising 
lattice, other solutions for regular two-dimensional lat- 
tices, such as triangularis honeycomb[(|, KagomeQ 
lattice and others were explored in several works and his 
importance in statistical physics waked up to search for a 
group of completely solvable models. The problem con- 
cerning to the exact solution and the critical behavior of 
the two-dimensional models was the matter of the Fan 
and Wu[l|-[§|. In those works the free fermion (FF) con- 
dition and the free fermion approximation was studied 
with great details and the relations of the Boltzmann 
weights for obtaining exact solvable models was estab- 
lished. In many situations when the FF condition is not 
satisfied completely it is possible to find with a good ap- 
proximation for particular values of the parameters of the 
model. This is the case investigated by Kun-Fa[l(j where 
the critical coupling of mixed Ising spin-1/2 with the ar- 
bitrary Ising spin-S was studied using the free fermion 
approximation. 

Since that many theoretical investigation was devel- 



oped, such as Ising-Heisenberg Kagome lattice[ll|, 12j. 
quantum square-Kagome antiferromagnetic lattice 13] . 
doubly decorated Ising-Heisenberg model[14j, the mixed- 
spin Ising model on a decorated square lattice with two 
different kinds of decorating spins on horizontal and ver- 
tical bonds [IB]- Other exactly solvable Ising model lattice 
known as square-hexagon (4-6) was considered by Lin an 
Yangfl6|. On the other hand a different 4-6 lattice as 
a special case of the 4-8 lattice was studied by Oitmaa 
and KeppertfljJ where the solution for the Ising spin-1/2 
case was found. It is remarkable to point out that the 
free fermion condition for the Boltzmann weights is sat- 
isfied identically in these models and the exact critical 
point can be performed, thus the models falls within the 
standard Ising universality class. 

Several real systems motivate to investigate in this 
kinds of lattice, such as the recently discovered 
two-dimensional magnetic materials Cu<jX^(cpa)^.xll-iO 
(cpa=2-carboxypentonic acid; X=F, CI, Br) where the 
Cu spins stands on the triangular Kagome lattice (lj| 
with Heisenberg interaction type. Liquid crystals net- 
works composed by pentagonal, square and triangu- 
lar cylinders Other recent investigation about the 
crystal structure of solvated [Zn(tpt)2/3(SiF,)(H20)2- 
(MeOH)] [tpt = 2,4,6-tris(4-pyridyl)-l,3,5-triazine] net- 
works with the (10,3)-a topology [20|. 



In this report we will discuss a similar model presented in ref. [17|, where we include local Heisenberg interaction. 
In FigQJi we display schematically two dimensional XXZ-Ising on 4-6 lattice, with a single line we represent the Ising 
lattice, whereas with double line we indicate the XXZ interaction. In fig. [TL we represent each decorated cell of the 
lattice, therefore the Hamiltonian of XXZ-Ising lattice will become as 



H({a},{S a }) = Y,JS i a i 
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the first summation runs over all sites involving the Pauli summation runs over nearest neighbor of whole lattice 
operator a with two possible values ±1, while the second containing the S a spin-1 /2 operators with a = {x, y, z}. 
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Figure 1: Schematic representation of two dimensional XXZ-Ising on 4-6 lattice. In (a) we represent as double line the XXZ 
interaction whereas with single line we represent the Ising interaction. In (b) we illustrate each decorated cell displayed in (a). 
The transformation of unitary cell is represented in (c), only as a function of the Ising interaction with nearest, next-nearest 
and quartic interaction of effective Ising lattice. 



The phase diagram at zero temperature for the two 
dimensional XXZ-Ising on 4-6 lattice is analyzed, com- 
puting the ground state energy. We verify that there are 
five different states, as depicted in figEl We can classify 
the states by ferromagnetic and antiferromagnetic in the 
Ising interaction part. Therefore let us start writing the 
ground state energy when the Ising interaction is ferro- 
magnetic, 



\FM) = 



++++ 
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|afi; 
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by the large ± (inner squared signals) we represent the 
Heisenberg interaction particles, whereas by the corner 
small signals ± we indicate the Ising interaction par- 
ticles. By R we represent the rotation operator acting 
only in Heisenberg interaction particles, each rotation is 
performed in ^ , around of the axis perpendicular to the 
plane of lattice. 

The ferromagnetic state (\FM)) given in eq.Q for 
J z < 0, is limited by |J 2 < J < 0, this state is de- 
picted as brown region in fig El There is also an anti- 
ferromagnetic state (|AFi)) represented by eq. ([3]) for 
J z < 0, which is restricted in the interval < J < — | J z , 
and is displayed in figEI as gray region. It is worth to 
notice that the state behaves as ferromagnetic for both 
in Ising and Heisenberg interaction particles but with 
opposite oriented spin. Under same condition we also 
have two types of ferrimagnetic states (FI) which are ex- 
plicitly given by eq.Q, for J z < 0, the \FI + ) is the 
ferrimagnetic state with magnetization 3/4, limited by 
f Jz < J < 0.43067503J Z , whereas the \FI~) corre- 
sponds to the ferrimagnetic state with magnetization 1/4 
restricted by -0.43067503 J z < J < -f J x . Those re- 
gions are illustrated in figEI as orange and cyan region 
respectively. These ferrimagnetic states are invariant un- 
der whole exchange of spin orientation. 



The other possible situation is when Ising interaction 
antiferromagnetic is considered. In this situation we only 
have one antiferromagnetic state \ AF 2 ), we also could call 
this state as frustrated statejH)]. This ground energy is 
obtained after diagonalized a 3 x 3 matrix, which is fall 
into a cubic equation and their lowest solution is given 
by — | (Pi cos(0i) — J z ), with 



p =v^FT9j2, 



(5) 
(6) 



in eq.((6]), 4>j (with j = 0, 1,2) is related to the real root of 
the cubic equation. The corresponding eigenvector state 
read as 



\AF 2 ) =hj^ R'' 



r=0 



(1 + 6 2 R) 



-++ 



and the coefficients of lO are given by 

h =^j(2P 1 cos(4> 1 ) + J z +3J) : 
b 2 =jjh (Pi Jcos(<£i) - J z ) - 1. 



(7) 



(8) 
(9) 



This antiferromagnetic state is present for arbitrary 
values of J when J z > 0, whereas for J z < this states is 
limited by I J| > -0.43067503 J z , as illustrated in fig. U 
We remark that in this case we have antiferromagnetic 
interaction for both Ising and XXZ interactions. 

It is worth to comment that this model has multicriti- 
cal point at zero temperature, where five states converge 
at J = J z = 0. 

To study the thermodynamics of the model we write 
the partition function of decorated XXZ-Ising lattice, 
given by Hamiltonian H]), 

Z(0)= J2 tr {s « } (e-^(M,{^})). (10) 

After taking the trace over operators {S a } we trans- 
form the decorated XXZ-Ising model into an effective 
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Figure 2: The phase diagram at zero temperature, as a func- 
tion of two parameters, J z and J. 



Ising model, with next nearest and quartic interactions 
parameter, whose effective Hamiltonian could be ex- 
pressed in general by 



H({a}) = K £ W+ L Y1 

ClCr2<73Cr4, 

<i,j> all 

square 

(11) 



with K being the nearest neighbor interaction, whereas 
L is the next nearest neighbor interaction parameter and 
M being the quartic interaction parameter. This trans- 
formation is also represented schematically in figHJ:. This 
effective Ising model is the so called 'Union Jack' lattice, 
which is an exactly solvable model[22|. 



Therefore the corresponding partition function of ef- 
fective Ising lattice is given by 



M=±i 



(12) 
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Figure 3: The eight spin vertex configurations. Reversal of 
all spins corresponds to the same vertex 

Using the Boltzmann weight given in figiU we are able 
to transform the Hamiltonian ([I]) into l|lip where their 
parameters are related by the following expression, 



/ =(wiw 2 wlwi) s , 
8/3 \w 2 
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Performing some algebraic manipulation we write the 
associated Boltzmann weights of the Hamiltonian[TJ This 
large results are written using some extra notations just 
to express in a compact form, 



=3 + c 4/3J ; + 8 cosh 3 (2/3J) + 2e" 4/3J ° cosh(4/3J) + 2q- 213 ^ cosh (2/3 y/ J 2 + 8 J 2 ) , 



w 2 =3 + 2c- 4 ' 3 ^ +4cosh(2/3J) +4008^2^/3^ +E e_S/3(PlC ° s(0jKJi) ) 

3=0 

w 3 =1 + c 4/3J - + 2e- 4/3Jz + 8cosh(2V2/3J) cosh(2/3J) + 2c^ J '(e 0P2 cosh(/3A+) + e~ pFii cosh(/3A_)), 

2 

w 5 =2 + cosh (4 cos (0 3 -) - 1)) + 2e 2/ ^ ( cosh(/3B + ) + cosh(/3B_)) + 2 cosh(2/3 J)(l + e~ 4/3 ^), 
j=o 



(17) 
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with Pi and <f>j already was defined in eqs.© and ([6]) 
respectively, whereas the other amounts are defined as 
follow 



2£j 

3 • 



p 2 =VJz + 10J 2 + 2Jt/4J2 + 25J 2 , 



At =^(3 ± ^)(J| + 4 J 2 ) + 12 J 2 - Pi 



B± =2d.P z + 5J 2 ± 2J^6J 2 + Jf . 



(21) 
(22) 
(23) 
(24) 



The other Boltzmann weights can be obtain using the 
symmetry rotation, thus we have the following identities 
for the model considered, 



W3 = Wi, and tog = we = wr = w$. 



(25) 



The two dimensional XXZ-Ising on 4-6 lattice model in 
general has no exact solution, but it could be possible to 
find some particular solutions imposing the exact solvable 
condition. Therefore it is possible to find one branch 
region where the model have an exact solution. 
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The first branch of a possible exact solution could be 
when the Boltzmann weights satisfy the so called sym- 
metric eight-vertex model (SEVM) condition, where we 
must have the following relations, 

Wi = W2, W3 = W4, W5 = wq and wj = wg. (26) 

Our model satisfy all these relations given by lf26"l) , un- 
less the first one. 




Figure 4: In (a) is displayed the exactly solvable condition 
for SEVM (wi = W2). (b) The FF is imposed and there is 
no exact solution. But there is a valley where lS./w'^ nax <C 1, 
particularly we show the region for A/iu^ a3; > 10 -3 as gray 
region. 

Imposing the first relation of eq. |26|) we have one pos- 
sible solution. In figQJa) we display the exactly solvable 
SEVM condition, as function of the parameters J and 
J z in units of /?. Therefore we show that for one free 
parameter the Hamiltonian (TJ could be solved exactly, 
spite this transcendental equation evolves complicated re- 
lation of J z and J, we are not able to invert one of them 
as a function of the other one explicitly, but even so we 
can invert numerically. In the limit for large values of J z 
and J, we have the asymptotic limit where the relation 
becomes approximately by J w ±2.713579J 2 . We also 
have a trivial solution when J = 0, this corresponds just 
to a set of non-interacting squared Ising model. 

The second candidate for the exact solution is the so 
called free fermion (FF) condition, when the following 
relation 

A = wiu'2 + W3W4 — W5W6 — (27) 

must satisfy the condition A = 0. 

Imposing the FF condition unfortunately we cannot 
find a solution for this model, unless for a trivial con- 
dition when J = 0. In this situation we always have a 
positive amount of A/w^ox- Then we can note that, if 
we display for small A / m^ax the relation of J z and J in 
units of (3, there are a wide valley where satisfy the con- 
dition A/to^ aa , <§C 1, particularly we display in fig.{4]), a 
valley for A/w^ a2 . < 10 -3 , whereas the gray region cor- 
respond to the condition of A/w^ax > 10 -3 . This means 



that we can approximate to the FF condition and solve 
this model with good approximation in all this region. 

It is also possible to discuss the critical behavior, even 
when exactly solvable condition is not satisfied. For the 
first branch solution (SEVM), its critical condition must 
satisfy the following relation 

W\ + W3 + u>5 + W7 = 2max(wi,W3,W5, wj). (28) 

In figHlwe display the critical points region as a function 
of the parameters J and J z in units of p, and we represent 
by a solid blue line. The convergence for this case is 
satisfied in all critical points \lSJ\/w^ nax < 1, with |A'| = 
\wi - w 2 \ and w max = max{wi,w 2 }. 

The second branch critical points region is when we 
impose the FF condition 

k'i + 102 + w 3 + w 4 = 2m&x(w 1: W2,w 3 , u> 4 ). (29) 

In figlUwe display the critical points region as a function 
of the parameters J and J z in units of (3, the case when 
the Boltzmann weight Wi is taken as the maximum value. 
The red solid lines indicates the region where the FF 
approximation is valid (\A\/wl < 1), and the doted red 
line indicates the region where \A\/wf > 1. The black 
solid line display the critical condition region when the w 2 
is the largest one, with restriction A/w| < 1, while with 
doted black line we represent the critical region when 
A/w| > 1. 




Figure 5: The critical points region under FF condition: red 
line corresponds the condition of wl > w2, whereas by black 
line we represent the condition when wl < w2. Doted line 
corresponds the region where A/w'^ nax > 1. On the other 
hand when SEVM condition is imposed, the critical region 
becomes the curve given by blue line. 

In this report we discuss some particular solution 
of two-dimensional XXZ-Ising model on square-hexagon 
lattice, where the decoration is a square with XXZ in- 
teraction and the interaction terms of the lattice is given 
by Ising type coupling. We discuss the phase diagram 
at zero temperature displaying five different phases. To 
study their thermodynamic, initially we used two param- 
eters but due to imposing the exact solvable condition we 
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constrain this two parameters, then we obtain a 2D XXZ- 
Ising on 4-6 lattice with one free parameter, under SEVM 
condition. Under FF condition, we display a wide valley 
where the model could be considered approximately as 
satisfying the FF condition. It is also possible to discuss 
the critical condition even when the exact result condi- 
tion is not satisfied. 
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